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Abstract. We attempt a unified treatment of heavy quarkonia and heavy- 
quark diffusion in the Quark-Gluon Plasma. Our approach is based on finite- 
temperature T-matrices with interaction potentials estimated from the heavy- 



■ quark internal energy computed in thermal lattice QCD (IQCD). In the charmo- 



nium sector S-wave bound states {J/ip, Vc) survive up to temperatures of ~2rc, 

■ not inconsistent with constraints from euclidean correlation functions in IQCD. 
00 ! In the open-heavy flavor sector, the T-matrix interaction reduces heavy-quark 

■ diffusion substantially, leading to fair agreement with single-electron spectra at 
^ I RHIC and suggestive for a small viscosity-to-entropy ratio close to Tc- 
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1. Introduction 

The properties of the strong-interaction matter produced in Au-Au collisions at 
RHIC are under intense debate. The empirical evidence collected thus far suggests 
a liquid-like medium which (i) quickly establishes, and then maintains, local ther- 
mal equilibrium for low transverse- momentum particles (pt<2 GeV), (ii) is largely 
opaque to high pt>5 GeV particles, (iii) exhibits constituent quark-scaling proper- 
ties in hadron production at intermediate pt~2-5 GeV. A key question is if and how 
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Fig. 1. Color-singlet heavy-quark internal energies at finite temperature (with the 
asymptotic value, Uoo=U{r oo,T) subtracted) extracted from 3-flavor [i8j (left 
panel) [14j and quenched IQCD [[9] (right panel) [[TO]. 

these phenomena are related [[TJ [2] and what the relevant interactions are (e.g., ra- 
diative or elastic, perturbative or nonperturbative, gluon- or quark-driven, etc.). In 
this article we approach this problem via the heavy-quark (HQ) sector. For heavy 
quarkonia, the hope has emerged that, using potential models, one can quantita- 
tively analyze in-medium quarkonium properties based on input from lattice QCD 
(IQCD), which is discussed in Sec. [21 The question arises whether similar interac- 
tions could be relevant for low-momentum properties of individual heavy quarks in 
the Quark-Gluon Plasma (QGP), i.e., for HQ diffusion; this is studied in Sec. [3l 
including applications to observables. We conclude in Sec. S) 



2. Heavy Quarkonia 

Our starting point is a two-body scattering equation for the in-medium T-matrix 
of a heavy quark (Q) and antiquark (Q) [131 H], 

T{E, q, q') = V{q, q') + J d'k V{q, k) G{E, k) T{E, fc, q') . (1) 

It represents a "ladder" resummation of a suitably defined HQ potential, V , and is 
widely used in the nuclear many-body problem as well as in the study of electromag- 
netic plasmas [[S]; G denotes the intermediate Q-Q propagator including selfenergy 
insertions due to interactions with the surrounding medium particles. The main 
approximations are neglecting (i) virtual Q-Q excitations, and (ii) retardation ef- 
fects in the interaction; both should be reasonable for large quark masses, mg. 
The T-matrix provides a uniform treatment of bound and scattering states which is 
particularly important for situations involving dissolving bound states, such as for 
atoms in plasma physics or quarkonia in the QGP. The precise relation between the 
HQ potential and free energy, F = U — TS, computed in thermal IQCD, is an open 
problem. We here adopt the internal energy, U, as the potential, but other choices 
are possible [ El [7] . Figure [T] shows that different extractions of U still imply up to 



Heavy Flavor in the sQGP 



3 




0.0 0.1 0.2 0.3 0.4 0.5 0.0 0.1 0.2 0.3 0.4 0.5 



T [fm] Z [fm] 

Fig. 2. Euclidean correlators for 77c (left) and Xci (right) in the QGP computed 
within the T-niatrix approach [ |3] using the HQ internal energy extracted from 
Nf—3 IQCD [|8]. The correlators are normalized to a reconstructed one based 
on the vacuum T-matrix. The increase in the Xc correlator is largely driven by 
zero-mode contributions [tllj implemented in quasi-particle approximation. 



~40% uncertainty. In either case, the increased color-screening leads to a reduction 
in charmonium binding with increasing temperature. However, when including an 
in-medium reduction of the c-quark mass governed by the asymptotic value of the 
internal energy, m*=m°-|-t/oo/2, the bound-state mass is almost stable. Together 
with the lowering of the cc threshold, -Ethr=2m*, and nonperturbative rescattering 
strength generated by the T-matrix [ [4] , the pertinent 77c correlator is rather stable 
(cf. left panel of Fig.©, similar to the findings in IQCD [[T11|T3]. The underlying 
spectral functions show that the bound state "melts" slightly above ~2Tc. Smaller 
melting temperatures are found in Refs. [[SI Ej based on different input potentials, 
but the IQCD correlators can be well reproduced, due to the above-mentioned in- 
terplay between binding and threshold effects. An independent determination of 
the in-medium HQ mass, an improved definition of the potential, a quantitative 
implementation of finite width effects, as well as a coupled channel treatment to 
account for gluonic excitations, may be required to make further progress. 

3. Heavy-Flavor Transport and Observables 

In Ref. [ the exchange of effective meson resonances in scattering of heavy 
quarks in the QGP has been introduced and found to reduce charm- and bottom- 
quark thermalization times by a factor of ~3 compared to elastic perturbative QCD 
(pQCD) scattering. When implemented into relativistic Langevin simulations for 
Au-Au collisions at RHIC [ 15 , the predictions for semileptonic electron (e^) spectra 
from HQ decays turned out to be in fair agreement with experiment [ 1161 I17j . 
In Ref. [ [18] , nonperturbative HQ interactions in the QGP have been evaluated 
by employing the same T-matrix approach as discussed in the previous section, 
including all color channels (1, 3, 6, 8) as well as S- and P-waves in heavy-light 
quark scattering. The pertinent T-matrices exhibit ( "pre-hadronic" ) resonance- 
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Fig. 3. Left panel: Imaginary part of tlie T-matrix for c-quark scattering of light 
quarks and antiquarks in the QGP Right panel: thermal relaxation rates for c- 
quarks following from the nonperturbative heavy-light T matrices (upper curves at 
p=0) and from LO pQCD scattering off quarks and gluons with as=QA. 



like D-meson and diquark correlations up to temperatures of ~1.5 and 1.2 T^, 
respectively, cf. left panel of Fig. [3] (using the potential in the right panel of Fig.[l] in 
the bottom sector, the dissolution temperatures are slightly larger). The repulsive 
sextet and octet channels, as well as P-waves, are suppressed. The dissolution 
of the resonances leads to a decrease of the thermalization rate with increasing 
temperature, opposite to the standard behavior as found, e.g., in pQCD, cf. right 
panel of Fig. [H 

The heavy-light quark T matrices (supplemented by pQCD interactions with 
gluons) have been used to compute HQ diffusion in a Fokker-Planck approach and 
applied to Au-Au collisions at RHIC utilizing relativistic Langevin simulations in 
an expanding fireball. The resulting nuclear modification factor and elliptic flow 
are quite comparable to the effective resonance model [[T31 [H] (see Fig. 01), with 
slightly smaller (larger) effects for c {h) quarks (which is qualitatively similar to 
the dissociation model of Ref. [HH])- Uncertainties due to different extractions of 
the IQCD internal energy [[lOl [205 amount to ^^30%. The HQ spectra have been 
hadronized into D and B mesons in a combined coalescence/fragmentation frame- 
work [IJH [15], with subsequent semileptonic (3-body) decays into single electrons. 
The relative weight of charm and bottom contributions is estimated from d-Au 
data [112], crossing at about Pj~5 GeV, (cf. also Ref. [H^)- Note that quark coales- 
cence at Tc naturally emerges from the "pre-hadronic" resonant correlations in the 
T-matrix [HI]- The electron spectra in Au-Au collisions [TS], which do not involve 
adjustable parameters, compare well to the most recent RHIC data [[THl [T71 [53], see 
left panel of Fig.[5l The HQ diffusion coefficient at p—Q can be used for a schematic 
estimate of the viscosity-to-entropy ratio {r}/s) in the QGP [fT]; evaluating their 
relation in both weak- and strong-coupling limits leads to the pink band in the right 
panel of Fig. [5] it rises with temperature and is indicative for a strongly coupled 
QGP close to Tc. 
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Fig. 4. HQ spectra in Au-Au collisions at RHIC utilizing relativistic Langevin 
simulations with nonperturbative HQ interactions [ I15|, I18j . 
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Fig. 5. Left: spectra in Au-Au collisions at RHIC following from the HQ spectra 
in Fig.[4]after hadronization and decay [[18|, compared to data [[TH [171 EH]. Right: 
estimates of rj/s [IT] for pQCD, resonance model and T-matrix approach. 
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4. Conclusions 

An in-medium T-matrix approach (utilizing potentials estimated from thermal lat- 
tice QCD) has been applied to evaluate HQ interactions in the QGP. For charmonia, 
S'-wave ground state can survive up to ~2Tc, roughly in line with IQCD correla- 
tors. For open heavy flavor, resonance-like correlations induce a small HQ diffusion 
coefficient which allows to describe data and suggests a small ry/s ratio. Several 
open problems remain, e.g., a proper definition of the potential, corrections to the 
T-matrix approach including radiative ones, and in-medium mass and width effects. 
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